Treatment of sound on quantum computers 
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We study numerically how a sound signal stored in a quantum computer can be recognized and 
restored with a minimal number of measurements in presence of random quantum gate errors. A 
method developed uses elements of MP3 sound compression and allows to recover human speech 
and sound of complex quantum wavefunctions. 
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In the last decade the rapid technological progress 
made possible the treatment of large amounts of informa- 
tion and their transmission over large distances. In spite 
of this the transmission of digital audio signals required 
the development of specific compression methods in order 
to achieve real time audio communication. A well known 
example of audio compression is the Mpeg Audio Layer 3 
(MPS) which allows to reduce the signal size by an order 
of magnitude without noticeable distortion [l]|. It essen- 
tially uses the Fast Fourier Transform (FFT) in order to 
have rapid access to the signal spectrum, whose analysis 
allows to reach significant compression rates. Such meth- 
ods find every day applications in the Internet telephone 
communication and teleconferences. 

Recent developments in quantum information make 
possible a new type of computation and communication 
using the quantum nature of the signal (see e.g. In 
quantum computation theory it was shown that certain 
quantum algorithms can be exponentially more efficient 
than any known classical counterpart. For instance the 
Shor algorithm enables to factorize large integers in a 
time polynomial in the number of bits whereas all known 
classical algorithms are exponential 0. This algorithm 
relies on the Quantum Fourier Transform (QFT) which 
is exponentially faster than FFT [3, 3 ■ Simple quantum 
algorithms have been realized experimentally with few 
qubit quantum computers based on Nuclear Magnetic 
Resonance (NMR) and ion traps 0,0, Quantum com- 
munications also attracted a great deal of attention since 
they allow to realize secret data transmission. Currently, 
such a transmission has been achieved over distances of 
up to a few tens kilometers 0. 

These prospects rise timely the question of treatment 
of audio signals on quantum computers. Classical audio 
analysis methods cannot be directly applied to quantum 
signals and it is important to adapt them to the new en- 
vironment of quantum computation. Furthermore while 
digital signal treatment is faultless, quantum computa- 
tion contains phase and amplitude errors which can affect 
the quality of sounds encoded on a quantum computer. 
In addition the extraction of quantum information re- 
lies on quantum measurements which bring new elements 



that must be taken into account in the treatment of quan- 
tum audio signals. Different type of sound signals are 
possible like human speech, music or pure quantum ob- 
jects like the Wigner function [|| which can be efficiently 
prepared on quantum computers 0, ^3 ■ 

For the standard audio sampling rate of AAkHz a quan- 
tum computer with 20 qubits (two-level quantum sys- 
tems, see 0) can store a mono audio signal of 23 sec- 
onds. A quantum computer with 50 qubits may store an 
amount of information exceeding all modern supercom- 
puter capacities (1000 years of sound). Thus the devel- 
opment of readout methods, which in presence of imper- 
fections can recognize and restore the sound signal via a 
minimal number of quantum measurements, becomes of 
primary importance. 

To study this problem we choose the following sound- 
track pronounced by HAL in the Kubrick movie "2001: 
a space odyssey": "Good afternoon, gentlemen. I am a 
HAL 9000 computer. I became operational at the H.A.L. 
lab in Urbana, Illinois on the 12th of January" [lll |. The 
duration of this recording is 26 seconds and at a sampling 
rate / = 8kHz it can be encoded in the wavefunction of a 
quantum computer with nq = 18 qubits (the HAL speech 
is thus zero padded to last 32 seconds). This rate gives 
good sound quality and is more appropriate for our nu- 
merical studies. Digital audio signals can be represented 
by a sequence of samples with values s„ in the interval 
(—1, 1) so that the n-th sample gives the sound at time 
t ~ n/ f. This signal can be encoded on a quantum com- 
puter by the following wavefunction = ^J2n^n\n >, 
where A is normalization constant. The state \n > repre- 
sents the multiqubit eigenstate |ai...ai...a„ > where 
is or 1 corresponding to the lower or upper qubit state, 
the sequence of gives the binary representation of n. 

Using numerical simulations we test various ap- 
proaches to the readout problem of the above signal en- 
coded in the wavefunction of a quantum computer. Di- 
rect measurements of the wavefunction do not allow to 
keep track of the sign of s„ and many measurements 
are required to determine the amplitude |s„| with good 
accuracy. Another strategy is to use the analogy with 
the MP3 coding. With this aim we divide the sound 
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(color online) Sound signal spectrum as a function of 



frame number x = n/An. On all panels the black curves show 
the spectrum of the original signal Sn- On the left top panel 
the green/gray curve represents the spectrum of s„ obtained 
with Ai = 5 measurements per frame in time domain. On 
the right top panel it shows the spectrum of s'„ for the same 
number of measurements M performed after QFT. Bottom 
panels show these spectra on a smaller scale. 



into consecutive frames of fixed size An = 2"^ where 
rif can be viewed as the number of qubits required to 
store one frame. We choose these qubits to be the nj 
least significant qubits in the binary representation of 
n = (ai...a„^_„j,+i...a„^). Then we perform QFT on 
these n / qubits that corresponds to applying FFT to all 
the 2'^''~^f frames of the signal in parallel. This requires 
nfijif + l)/2 quantum gates contrary to 0{nf2"f) clas- 
sical operations for FFT. After this transformation the 
wavefunction represents the instantaneous spectrum of 
the sound signal evolving in time from one frame to an- 



other. This way the most significant Uq 



qubits store 



the frame number k while the least significant Uf qubits 
give the frequency harmonic number j. Hence, after QFT 
the wave function has the form ip — "^f^ j Sk.j\k, j > 
where Skj is the complex amplitude of the j-th harmonic 
in the k-th frame. The measurements in this represen- 
tation gives the amplitudes \Sk.j\ while phase informa- 
tion is lost. However for sound the main information is 
stored in the spectrum amplitudes and the ear can re- 
cover the original speech even if the phases are all set to 
zero. Thus the recovered signal is obtained by the inverse 
classical FFT and is given by s',, = J2j 15"^^ je^'^^J™/^" 
with n = k An + m (to listen sound we use i?e(sjj)). For 
time domain measurements of s„ the recovered signal is 
Sn = \sn\- These expressions for and s„ hold for an in- 
finite number of measurements. In reality it is important 
to approximate them accurately with a minimal number 
of measurements. 

For our soundtrack we found that the optimal frame 
size is An = 2^(n/ = 9) and we perform M measure- 
ments per frame. In Fig. 1 we compare the spectrum 
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FIG. 2: (color) Coarse grained diagram for the signal of Fig. 
1 obtained by measurements of qubits 1 to 4 and 10 to 13 
(see text). Top left: original signal; bottom left: distribution 
obtained after QFT with the same number of measurements 
as in Fig. 1; bottom right: same as bottom left with the noise 
in the quantum gates (e = 0.05); top right: diagram obtained 
from measurements in time domain (with the same number 
of measurements as in bottom right). The color represents 
amplitude of the spectrum: blue for zero and red for maximal 
values. The horizontal/ vertical axis corresponds to coarse 
grained frequency /time. 



of s'„ and s„ with the original signal spectrum. Here 
only M = 5 measurements per frame are performed and 
the results clearly show that the quality of the restored 
sound is significantly higher for the spectrum domain 
measurements. Examples of restored and original sounds 
are available at The HAL speech is recognizable 

from Af = 5 for spectrum domain measurements while 
it is distorted beyond recognition for direct time domain 
measurements even for AI — 100. 

Fig. 1 shows the global structure of the signal spec- 
trum. To make comparison more quantitative and visual 
we show coarse grained color diagrams of the spectrum. 
The coarse graining is obtained by measuring only cer- 
tain qubits corresponding for example to aia2a3a4 and 
aioaiiai2ai3. For this gives a coarse grained spec- 
trum \Skj \ in 2* X 2^ cells shown in Fig. 2. The same 
total number of measurements as in Fig. 1 allows to 
reproduce the original coarse grained diagram with good 
accuracy. Even if QFT is performed with noisy gates (the 
angle in the unitary rotations fluctuates with an ampli- 
tude ETT = O.OStt) the spectrum diagram remains stable 
and is reproduced with good accuracy (see Fig. 2). At 
the same time the coarse graining in the time domain for 
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FIG. 3: Fidelity / for sound signals s'^ (open circles) and s„ 
(full circles) as a function of number of measurements M. For 
sjj fidelity is shown for various amplitudes of noise in the QFT 
gates with e = (o), 0.05 (diamonds), 0.1 (triangles), 0.3 (+), 
1 (squares). Inset shows data at small M scale. 
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FIG. 4: Dependence of fidelity / for the sound signal on the 
noise amplitude e for the number of measurements M — 1000 
(o), 100 (squares), 20 (triangles), 5 (diamonds). Inset shows 
a data fit 1 - / = 0.36e^n^ (full curve) at smaU e, M = 100. 



the signal Sn shown in Fig. 1 gives the diagram which is 
very different from the original. 

The global quality of the recovered signal s„ 
(or s^) obtained via a finite number of measure- 
ments is convenient to characterize by the fidelity 
defined as / = \J2n^n^*Sn\M,e)\/R with R = 

(E„|sn^l'E„ki^(M,e)n^/'. Here, sV{M,e) is the 
signal obtained in a way described above with M mea- 
surements per frame in time domain (s„(M)) or in fre- 
quency domain after QFT with noisy gates {s'„{M,e)). 
The dependence of / on M is shown in Fig. 3. For 
large M it approaches to unity for both signals s„ and 
at e = 0. However, for a small number of measure- 
ments (5 < M < 50) the fidelity is significantly higher 
for measurements performed in the frequency domain af- 
ter QFT (Fig. 3 inset). The presence of noise in the 
quantum gates used in QFT for s'^ reduces the value of 
/ but for 5 < M < 50 and e < 0.1 this reduction is not 
significant. The drop of / becomes considerable only at 
relatively large amplitudes with e > 0.2 as it is shown 
in Fig. 4. The residual level of / at maximal e « 1 
is in agreement with the statistical estimate according 
to which /(e = 1) w ^JuifWi w 0.2, where in is the 
number of frequencies per frame for the original signal 
(ui K, 20 according to Fig. 2). At small e the drop of / 
is quadratic in e (1 — / ~ e^n^) since each gate trans- 
fers about of amount of probability from ideal com- 
putational state to all other states jl(J|. The obtained 
results show that the MP3-like strategy adapted to the 
quantum signals allows to recover human speech with a 
significantly smaller number of measurements with a re- 
duction factor of 10-20. 

Above we assumed that the sound signal is already en- 
coded in the wavefunction. For certain quantum objects 



such an encoding can be done efficiently. As an example 
we consider the wavefunction evolution described by the 
quantum sawtooth map 

where I — —id/dO, h = 1, k,T are dimensionless map 
parameter and ip is the value of ip after one map it- 
eration (we set h — 1). In the semiclassical limit 
k ^ 1, T <C 1 the chaos parameter of the model is 
K — kT — const. The efficient quantum algorithm for 
the simulation of this complex dynamics was developed 
and tested in [isl . The computation is done for the 
wavefunction ^ on a discrete grid with N = 2"' points 
with 9n ~ 2Tm/N, n = 1, ...,7V in ^-representation and 
I + iV/2 — 1, .., iV in momentum representation. Here, as 
before Uq is the number of qubits in a quantum computer 
and in ^-representation ij^ = X^n '4'iPn)\'n > is encoded in 
the register \n >~ |ai...ai...a„^ >. The transition be- 
tween n and 9 representations is done by QFT and one 
map iteration is computed in O(n^) quantum gates for an 
exponentially large vector of size 2"' To study the 
sound of quantum wavefunctions of map ^ we choose 
here a case with K = —0.5, T = 2tt/N and Ug = 14 
corresponding to a complex phase space structure. 

The signal encoded in the wavefunction ip{(^n) after t 
map iterations can be treated in a way similar to one 
used before for the HAL speech s„. The measurements 
in 0-basis give the signal s„ = |'!/'(^n)| which however 
requires a large number of them to suppress noise (also 
the phase is completely lost). Another method works as 
for signal: first QFT is performed on Uf = 5 less sig- 
nificant qubits giving tp = J2k j ^k,j\k,j > and then the 
measurements are done to determine the instantaneous 
spectrum amplitudes \Skj\ of ip{9n) (here k = 1, 2^ is 
the frame number, j = 1, 2^ is the index of frequency 
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FIG. 5; (color) Left: coarse grained diagram of S'®' for the 
sound signal s'n obtained from the quantum computation Q 
of the wavefunction after t — 100 map iterations, top panel 
shows exact distribution S'®' and bottom panel is for a num- 
ber of measurements as in Figs. 1,2. Bottom right: same as 
bottom left but with noise amplitude e — 0.05 in the quantum 
gates. Right top: the exact Husimi distribution h{l,9). The 
initial state is a momentum eigenstate at I — 100 — N/2. The 
color and axes are as in Fig. 2. 



harmonics and An = 2^). The sound of quantum wave- 
function is recovered via the inverse classical FFT giving 

signal defined before. Examples of restored sound are 
given at |0| and clearly show that the quality of MP3- 
like signal is much higher compared to s„ (we use 
sampling rate f — 1 kHz for this case with Uq = 14). 

A more detailed analysis of the quantum sound can 
be obtained from the coarse grained spectrum digram 
similar to the one in Fig. 2. The coarse graining of Skj 
is done by measuring 5 most significant and 5 less signif- 
icant qubits corresponding to aia2...a5 and aioaii...ai4 
that gives coarse grained distribution S^^^ in 32 x 32 cells. 
The diagrams of S^^^ obtained for infinite and finite num- 
ber of measurements are displayed in Fig. 5 (left top and 
bottom respectively) . The exact diagram shows an inter- 
esting structure which is recovered with a finite number 
of measurements. This structure remains robust against 
noise in the quantum gates used for computation of 100 
map iterations and final QFT (Fig. 5 right bottom). 

The origin of this structure becomes clear after its com- 



parison with the coarse grained Wigner function called 
the Husimi distribution h{9, 1) 0,E3| which is shown in 
Fig. 5 (right top) which is very close to (left top). 
Indeed, h{9J) is defined in the phase space (1,9) by 



where the gaussian smoothing function is G{1' — I) = 
(T/7r)i/4e-^(''-')'/7ViV [H m. The Husimi distri- 
bution is always positive and gives a direct compari- 
son between the classical phase space Liouville density 
distribution and a quantum wavefunction. In fact the 
coarse grained distribution 5*'^' is also given by equation 
^ where G{V) is replaced by a constant in the interval 
Al' = 2^ and zero outside that corresponds to the appli- 
cation of QFT to less significant qubits n/ = 5. Such a 
replacement modifies the values of coarse grained h{9, 1) 
but this modification remains small if A?' > 1 [3. As 
a result we may argue that the signal represents the 
quantum sound of coarse grained Wigner function. 

In conclusion, our results show that sound signals 
stored in a quantum memory can be reliably recognized 
and recovered on realistic quantum computers. The 
method proposed allows to obtain sound of quantum 
wavefunctions that can be useful for future quantum 
telecommunications . 
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contract No. DAAD19-01-1-0553. 
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